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PACS 83.60. Df “ Nonlinear viscoelasticity 

PACS 83.80.Hj - Suspensions, dispersions, pastes, slurries, colloid 
PACS 83.80. Fg - Granular solids 

Abstract - When pulling a probe particle in a many-particle system with fixed velocity, the 
probe’s effective friction, defined as average pulling force over its velocity, 7 e// := {F^x)/u, first 
keeps constant (linear response), then decreases (thinning) and finally increases (thickening). We 
propose a three-time-scales picture (TTSP) to unify thinning and thickening behaviour. The 
points of the TTSP are that there are three distinct time scales of bath particles: diffusion, 
damping, and single probe-bath (P-B) collision; the dominating time scales, which are controlled 
by the pulling velocity, determine the behaviour of the probe’s friction. We confirm the TTSP by 
Langevin dynamics simulation. Microscopically, we find that for computing the effective friction, 
Maxwellian distribution of bath particles’ velocities works in low Reynolds number (Re) but fails 
in high Re. It can be understood based on the microscopic mechanism of thickening obtained in 
the T — 0 limit. Based on the TTSP, we explain different thinning and thickening observations 
in some earlier literature. 
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Introduction. — Microrheology studies flow of com¬ 
plex fluids under micro-mechanical control [^[^. It pro¬ 
vides not only a novel method to understand materials’ 
viscoelasticity on the microscopic level but also a 
nice example of studying the response theory, which is a 
fundamental issue in statistical mechanics [7|-[l0]. While 
in passive microrheology, only linear response to thermal 
fluctuation is possible, in active microrheology (AM), non¬ 
linear response can also be realized by large pulling. A 
typical AM experiment is to pull a probe particle em¬ 
bedded in a complex fluid with fixed velocity and then 
measure the fluctuating force of it The probe’s fric¬ 
tion coefficient, defined as average force over its velocity 
7 e// ■“ {Fex)/u, shows intriguing nonlinear behaviour: it 
first keeps constant (linear response) in the small velocity 
regime, then starts to decrease (thinning) in the moderate 
velocity regime, and finally may increase (thickening) in 
the large velocity regime. Similar behaviour can occur in 
bulk shear of macrorheology 


11-13 as well. 


Linear response and thinning were observed in col¬ 
loidal systems both in experiments 
tions 
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body Smoluchowski equation in low density 18 , density 


^One can also pull the probe with constant force, then measure 
the fluctuating velocity. 



a unifying description of thinning and thickening is still 
absent. Why was only thinning observed in colloidal sys¬ 
tems, while only thickening was observed in static granular 
systems? In this letter, we address this issue by proposing 
a TTSP. 

Three-Time-Scales Picture. — A unifying picture 
of thinning and thickening is that three time scales of bath 
particles are involved (see fig. [^): 

• diffusion time scale: t^iff = R^/D, where D = 

is the diffusion coefficient with the solvent friction 70 , 
R is the characteristic length scale (for hard sphere 
systems, it should be the center distance of the probe- 
bath particles contacting with each other). The cor¬ 
responding diffusion velocity is Udiff = R/tdiff- 

• damping time scale: tdamp = nT-b/jo, where to;, is 
the mass of a bath particle. The damping velocity is 

^damp — R/td amp ■ 
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• collision time scale: tcoi = R/u, where u is the pulling 
velocity. It characterizes the mean-free time between 
first and second P-B collisions without damping. The 
collision velocity is Ucoi = R/tcoi = u. 


The dominating time scales are controlled by the pulling 
velocity u, which can be indicated by Peclet number Pe := 
u/udiff = and Reynolds number Re := ujudamp = 

Different dominating time scales lead to different 
behaviour of the increased friction Aye// = 7e// — 7o- 
In detail, (i) when the pulling velocity is small enough 
that Pe <C 1 and Re <C I, the diffusion dominates. Aye// 
arises from the diffusion of bath particles, which leads to a 
linear response regime, (ii) As the pulling velocity is much 
larger than the diffusion velocity but still much smaller 
than the damping velocity, i.e. Pe ^ 1 and Re 1, dif¬ 
fusion is unimportant, damping dominates. Aye// arises 
from the damping of bath particles, which leads to an¬ 
other linear response regime, (iii) As the pulling velocity 
is even larger than the damping velocity, i. e. Pe ^ 1 and 
Re 3> 1, inertia dominates, Aye// arises from single P-B 
collision, which leads to an increasing friction regime. 

The plateau value of the linear response regime in (i) 
should be larger than the value in (ii), because diffusion 
causes larger friction in (i) comparing to the one arising 
from the damping only in (ii). As a result, the crossover 
from (i) to (ii) causes thinning. And the crossover from (ii) 
to (iii) causes thickening. The turning points of thinning 
and thickening should be around Pe = 1 and Re = 1, 
respectively (see fig.[^. 


Model. — To demonstrate the TTSP, we consider the 
model of pulling a probe particle with fixed velocity em¬ 
bedded in a suspension of N identical bath particles in two 
dimensions. Because pulling with fixed force and pulling 
with fixed velocity behave similarly, both may show thin¬ 
ning and thickening]^ we choose the latter for simplicity. 
All particles are assumed to be smooth and elastic hard 
disks with the same radius rp. The dynamics of a bath 
particle (labelled i) and of the probe (labelled p) obey the 
Langevin equations (la I and respectively. 


T ^2 "f- Fiscal 

0 — T T Fp.col T Fex (It)) 


where nib is the mass of a bath particle; Vi is the ve¬ 
locity of the i-th bath particle, u is the fixed pulling ve¬ 
locity of the probe; yo is the friction coefficient (all par¬ 
ticles have the same value due to yo oc roy; rj is the 
solvent’s viscosity.); £,k {k = i or p) is a Gaussian ran¬ 
dom force satisfying the fluctuation-dissipation relation 
(^fc(t)^fc,(t')) = ‘2lokBTSk,k'S''’^S{t-t') {x,y} are 

the components of the random force); Fiscal (or Fp^coi) is 
the interaction force between particles; F^x is the external 

^The effective friction of pulling with fixed velocity in gener al is 
larger than the one of pulling with fixed force as pointed out in ||18] 
and further analysed in [27| 



Fig. 1: Sketch of the TTSP of thinning and thickening: the 
effective friction ye// = F^x/uvs. the pulling velocity «. Three 
time scales of bath particles are involved: diffusion, damping 
and single P-B collision, which lead to three friction behaviour: 
a high plateau regime, a low plateau regime, and an increasing 
friction regime, respectively. The pulling velocity controls the 
dominating time scales. The crossovers cause thinning and 
thickening. The turning points of thinning and thickening are 
around Pe := ujudifs ~ ~ ^ u/udamp = 

= 1, respectively. 


pulling force on the probe only. According to the equa¬ 
tion of motion (EOM) (Ib), the probe’s increased friction 
^7e// := leff - 7o is 


^7e// — {Fp^col)/'^ 


( 2 ) 


where m, Fp ^oi are absolute values of the correspond¬ 
ing vectors, and the random force has been averaged out: 
(Vp) = 0 . 

Obviously, the P-B interaction directly leads to the in¬ 
creased friction, while the bath-bath (B-B) particle inter¬ 
action affects indirectly. We omit the B-B interaction in 
our model because 1) such interaction may not be neces¬ 
sary for thinning and thickening behaviour; 2) the omis¬ 
sion itself should be valid in the low density limit. In 
addition, we set the mass of the probe as much heavier 
than the mass of a bath particle: mp ^ mf,, so that in the 
coordinate of the probe, P-B collision just causes specu¬ 
lar reflection of the bath particles, but affects little the 
probe’s velocity. 

Now the system of pulling a probe with fixed velocity u 
is equivalent to the system of a flow with velocity —u of 
a suspension of N non-interacting bath particles passing 
a fixed disk with radius R = 2ro. The EOM of a bath 
particle (the index i is dropped) in the coordinate of the 
probe is 

m/d =-yo(u-I- m)-I- 4 (3a) 

with the reflecting boundary condition (RBC) 

V = V — {v ■ e„)e„ for |r| = R, (3b) 


where R = 2ro is the contact distance between the probe 
and a bath particle, and e„ is the unit normal vector along 
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Fig. 2: (a) Simulation results for the increased friction Aye// vs. the pulling velocity u for different solvent frictions and 

temperatnres. Different colors in (b) and (c) indicate different solvent frictions and temperatnres as these labelled in (a), (b) 

The rescaled increased friction Aye/z/yo vs. the Peclet nnmber Pe = u/udiff = The red arrow indicates the turning 

point of the thinning around Pe = 1. At Pe = 100, for different solvent frictions and temperatures, the corresponding value of 
the Re numbers are indicated, (c) The rescaled increased friction Aye///yo vs. the Reynolds number Re = ujudamp = 

The red arrow indicates the turning point of the thickening around Re = 1. At Re = 100, for different solvent frictions and 
temperatures, the corresponding value of the Pe numbers are indicated. On the plot of {yo = 100, ksT = 0}, four filled blue 
circles at Re = 0.01, 0.3, 1, 5 are drawn to compare with the corresponding streamlines in hg. 


the direction from the center of the probe to that of the 
bath particle colliding with it. Note that the P-B interac¬ 


tion term Fiscal in Eq.(la) is mapped into the RBC 
Being equivalent to its stochastic description Eq.([^, the 
probability description of a bath particle obeys the Eokker- 
Planck equation (EPE) 


dtP(r,v,t) = —v-drP+—dv- 
mb 


{v+u) + ^^^dj, 


mb 


P, (4a) 


which can be obtained by Kramers-Moyal expansion [28| 
of Eq. ([^. The corresponding RBC is 


P(r, V, t) = P{r, V - {v ■ e„)e„, t) for |r| = R. (4b) 

In principle, the steady state equation {dtP = 0) of the 
EPE (4a) can be solved with the RBC (|4b|). Then one can 


obtain the average collision force of N bath particles on 
the probe: 

{Fp,coi) = j dv j) dlNPst{r,v)v ■ (-e„)0[t; • (-e„)] 
2{-en)mbV ■ (-e„) 

( 5 ) 

where Pst denotes the steady distribution, dlNPst{r, v)v ■ 
(—e„)0[i; • (—e„)] is the density current of bath particles 
with velocity v passing through a small contact surface 
dlEn ( 0(a;) = 1 for X > O;0(x) = 0 for x < 0 ), and 
2(—e„)TOhU- (—e„) is the bath particles’ momentum trans¬ 
ferred to the probe due to single P-B collision. Inserting 
Eq.([^ into Eq.([^, one obtains the effective friction ye//- 
In practice, however, to analytically solve the EPE ([^ 
is difficult due to the RBC (4b). Our strategy is to numer¬ 
ically solve Eq. @ by simulation of the stochastic process 


Eq.([^, because of its equivalence to the EPE (|^ and sim¬ 
plicity. 

Stochastic Simulation. — To calculate the effective 
friction, the stochastic dynamics simulation according to 
Eq. (§ is performed. The discrete form of the Gaussian 
random force is ^ = y/2^ where (/i S 

X, y) is the standard Gaussian random number of the prob- 

ability distribution function as P{^q) = 2 )’ 

and h is the time step of the dynamics set to be h = ^ 
for different solvent frictions. The box size is set to be 
Lx X Ly = 20R x 201? with periodic boundary conditions, 
which is large enough to suppress finite size effects. The 
mass of bath particles and the P-B contact distance are 
set to be unit values: Wp = 1, i? = I. The density of bath 
particles is also rescaled to unit value no = I, since it is not 
a control parameter in our model due to the assumption 
of non-interacting of bath particles. The control parame¬ 
ters are the pulling velocity u, the solvent friction yo and 
the temperature ksT, which are applied to investigate the 
whole regime of different time scales. 

Initially, bath particles are homogeneously distributed 
in space with Maxwellian distributed velocities. Then 
the probe is pulled along the x direction with total run¬ 
ning time lOR/u, which ensures that the bath particles 
around the probe reach the steady state. After a tran¬ 
sient time, the steady average P-B collision force is com¬ 
puted by detecting the bath particles passing through the 
boundary: (Fp^coi) = dt2[mp{-en)v{t) ■ 

(—e„)]0[u(t) • (—e„)], which is the simulation realization 
of the collision force expressed in Eq. (§. The corre- 
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spending increased friction is obtained based on 

Eq. (H. 

Result. — Fig. (a) shows the simulation result 
of the increased friction Ajefj versus the pulling ve¬ 
locity u for different solvent frictions and temperatures, 
{70 = 1000 , ksT = 1000 , 100 } and {70 = 100 , ksT = 
100,10,0}. All plots, except for {70 = 100, ksT — 0}, ex¬ 
hibit linear response, thinning and thickening as expected 
by the TTSP. For the exception, only linear response and 
thickening occur, because no diffusion but only damping 
and collision time scales are involved. 

Fig. I (b) 0 shows that the rescaled increased friction 
Ajeff/jo versus Peclet number, Pe = u/udiff = 

In the small Pe regime Pe < 1, the diffusion time scale 
dominates, all plots coincide with each other in a plateau 
value. With increasing Pe, diffusion becomes less impor¬ 
tant, all plots start to decrease around Pe = 1, which 
agrees with the TTSP. Between Pe = 100 and Pe = 1000, 
for {70 = 1000, ksT = 100}, the brown line, clearly there 
is a second plateau lower than the first one, being consis¬ 
tent with the TTSP. In addition, the length of the thinning 
regime varies for different data sets|^ because for the same 
Pe, the Re numbers can also be different. At Pe = 100, 
for {70 = 1000, ksT = 100}, Re = 0.01, bath particles are 
still in the damping regime; for {70 = 100, ksT = 100}, 
Re = 1, bath particles are already in the inertia (thicken¬ 
ing) regime, which suppresses the thinning process. 

Fig. (c) shows that the rescaled increased friction 
A'yeff/lQ versus Reynolds number, Re = u/udamp = 
All plots start to converge around Re = 1, which 
agrees with the TTSP. In the small Re regime, for dif¬ 
ferent plots, at Re = 0.01, the frictions increase with the 
decreasing Pe as indicated in the figure, which supports 
the TTSP that the diffusion causing larger friction than 
the one in the damping only regime Pe —>■ 00. For Re > 1, 
i.e. the inertia regime, all plots coincide with each other 
and asymptotically tend to Aye// oc u, because the flux 
of bath particles passing through the P-B contact sur¬ 
face is j oc uqu with momentum transferring to the probe 
p oc mu, and Aye// = {Pcoi)/u = jp/u oc u. 

In summary, the friction behaviour of different domi¬ 
nating time scales and of the two turning points as shown 
in fig. all agree quite well with the TTSP. 

Microscopic picture. — The TTSP is indicated by 
Pe and Re. Microscopically, what happens in the different 
Pe and Re regimes? 

density distribution. It is convenient to compare the 
behaviour of bath particles in different Pe regimes by com¬ 
puting the pair distribution function g{r), which is the 
normalized number density of bath particles in the coor¬ 
dinate of the probe g{r) = n(r)/no = V J dvp{r,v) (V 

®Data set of {70 = 100, kgT = 0} is not included in fig[^(b), 
because u^i // = 0 , no diffusion is involved. 

^Tlie exception is that {70 = 1000, fcsT = 1000} and {70 = 
100, kgT = 10} coincide with each other in both figm (b) and (c), 
because for the same u, they have the same Pe and Re numbers. 


is the volume of bath particles, here for 2 d it is the area, 
no = N/V). Fig. 1^ shows the simulation result g{r) for 
different Pe numbers. For small Pe numbers Pe = 0.1,1, 
bath particles are both built up in front and left behind 
of the probe, i.e. the diffusion dominating regime. As Pe 
is quite large, Pe = 10,100, only a thin layer of bath par¬ 
ticles build in front but no particles left behind in a long 
tail region of the probe, which means that the diffusion 
is ignorable. The observation that diffusion dominates in 
the small Pe and is unimportant in large Pe, is consistent 
with the TTSP. 



Fig. 3: Pair distribution function g(r) of bath particles in the 
frame of the probe for different Peclet numbers. The solvent 
friction is 70 = 100 and the temperature is ksT = 100 . 


velocity distribution. Does the pair distribution func¬ 
tion contain enough information to calculate the effective 
friction? If the velocity of bath particles is Maxwellian 

distributed: feq(v) = (with thermal ve¬ 
locity Vth = s/ksT/m}y, then the total probability can be 
separated into p{r,v) = V~^g{r)feq{v), and the collision 
force in Eq. ([^ is reduced to 

{Fp^coi) = -noksT 0 dleng(r), ( 6 ) 

Jr^R 

which is identical to the one in ref. [^. If the velocity is 
delta distributed, f{v) = 5{y — {—vYf, the collision force 
in Eq. © is reduced to 

{Fp,coi) = -no2m& f dleng{r){u ■ e„)^0[M • e„] (7) 
Jr=R 

Inputting g{r) from the simulation into Eq.Q and Eq.Q, 
respectively, we obtain two increased effective frictions, 
the green line and black line, respectively, as indicated in 
fig-i Comparing them with the direct simulation result, 
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Fig. 4: Increased effective friction from different methods: 
the violet line is the direct simulation result, the green and the 
black lines are calculated by inputting g{r) from the simulation 
associated with Maxwellian distribution and delta distribution 
of velocity parts, respectively. 



the violet line (see fig. , one can conclude that for the 
calculation of the friction, the pair distribution function 
still works, but the proper velocity distributions should 
be input, according to different Re regimes: Maxwellian 
distribution in low Re and the delta distribution 6(v — 
{—u)) in high Re. 


streamline at T = 0. To further investigate the role 
of the Reynolds number, let us consider the T = 0 limit, 
where the diffusion time scale is ruled out, Pe —> oo. 
Eq.([3a|) is reduced to 


TOfcb = -"fov - 7oM 


( 8 ) 


with the RBC (3b). Interestingly, such simple dynamics 
provides a clear mechanism of thickening: the crossover 
from creep flow in the low Re to gas-like (inertial) flow in 
high Re, see fig. The black curves are the streamlines 
of the bath particles in the frame of the probe; red arrows 
are the velocity field. Before any collision, bath particles 
are moving with a constant velocity —it. Collision causes 
mirror-like reflection. The term —jqV reduces the veloc¬ 
ity, while —7oit accelerates it. A loose criteria of single¬ 
collision-only should be utdamp > R, be. Re = > 1- 

In the small Re limit, many P-B collisions occur and the 
bath particles tend to creep along the surface, see fig. 
Re = O.OI, which causes Fcoi oc u and Aye// oc i(°(the 
proof will be given somewhere), while in the large Re limit, 
the single collision causes Fcoi oc and Aye// oc it^. 

Based on the microscopic picture of Re (fig. [^, we can 
also understand why Maxwellian distribution works in low 
Re but fails in high Re. Let us consider the Pe 3> 1 
limit, where a bath particle moves with velocity —it rela¬ 
tive to the probe before any P-B collisionj^ Re determines 

® Indeed, before any P-B collision, the motion of a bath particle 
is determined by Pe only. It has nothing to do with Re. 


Fig. 5: Streamlines of bath particles for different Reynolds 
numbers at the T = 0 limit. Red arrows are the velocity field 
(rescaled by the pulling velocity for comparison). For small 
Reynolds number, the bath particles behave like creep flow 
around the contact surface. For high Reynolds number, they 
behave gas-like: single collision and flying away. The corre¬ 
sponding increased frictions are indicated with four filled blue 
circles in the gray line in fig. [^c). 


whether the solvent plays a role during P-B collisions. 1) If 
Re <IC I, damping dominates, the injecting velocity —it of 
the bath particle is quickly ’’erased” due to the damping 
and agitation processes by the solvent at the beginning 
of a few P-B collisions. In the following many times P- 
B collisions, the bath particle transfers the thermalized 
velocities to the probe. That’s why Maxwellian distribu¬ 
tion works in this limit, b) If Re ^ 1, inertia dominates, 
the P-B collision happens once only. The bath particle’s 
velocity transferring to the probe is exactly the injecting 
velocity —ii, which has nothing to do with the solvent. 
Thus, instead of the Maxwellian, the delta distribution 
6(^v — (—It)) works in this limit. 

Conclusion. — We propose a TTSP to unify thinning 
and thickening phenomena in active microrheology (see 
fig-0, and confirm it by a model of pulling with fixed ve¬ 
locity. The simulation result (fig. 0, which is equivalent 
to the solution of the FPE (0 in steady state, shows linear 
response, thinning and thickening. As far as we know, this 
is the first example demonstrating that both thinning and 
thickening can occur in non-interacting bath particles sys¬ 
tems (only P-B interaction is included), which indicates 
that the many body interaction is not necessary for thin¬ 
ning/thickening behaviour in the low density. Further- 
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more, as shown in fig. the results of the turning points 
of the thinning and thickening being around Pe = 1 and 
Re = 1, respectively, and the friction behaviour in differ¬ 
ent time scale regimes, all agree with the TTSP. 

Microscopically, the pair distribution function g{r) is 
obtained from the simulation as shown in figs|^ For the 
calculation of the friction, we find that with the input 
g{r) from the simulation, Maxwellian distribution works 
in low Re, but fails in high Re; while the delta distribution 
6(v — (—m)) works in high Re, but fails in low Re. In 
the T=0 limit {Pe —>■ oo), we obtain a clear microscopic 
picture of thickening for different Re regimes, see figj^ 
When Re ^ 1, damping dominates, the constant friction 
comes from creep flow, the bath particles collide with the 
probe and then creep around it; when Re 3> 1, inertial 
dominates, the increasing friction comes from the single 
P-B collision. Based on the picture of bath particles in 
different Re regimes, the validity/invalidity of Maxwellian 
distribution can also be understood. 

According to the TTSP, thinning arises from the 
crossover from diffusion to damping, and thickening arises 
from the crossover from damping to inertia. Note that 
diffusion was not involved in the experiments of pulling 
a single particle in static (T = 0) granular systems 
23 24 , that’s why thinning was not observed. For the 
same reason, it was not included in our earlier kinetic 
Thickening was not found in colloidal sys- 


model 

terns 


26 


14 16 18 19 , because they were limited to i?e <C 1 


regime, where inertia is unimportant. 

The TTSP should also be valid in the high density with 
dressed values of Pe and Re. B-B many body interaction 
increases the friction of a single bath particle, yg > yg 
(in the low density limit, yg is just the solvent friction 
yg). Based on the TTSP, the turning point of thinning 
Pe = 1 (X uyg should shift to a smaller pulling velocity 
value, and that of the thickening Re = 1 oc u/yg should 
shift to a larger value. 
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